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A double-exchange mechanism for the emergence of ferromagnetism in cubic uranium compounds is 
proposed on the basis of a j-j coupling scheme. The idea is orbital-dependent duality of 5/ electrons 
concerning itinerant Tg and localized iy states in the cubic structure. Since orbital degree of freedom is 
still active in the ferromagnetic phase, orbital-related quantum critical phenomenon is expected to appear. 
In fact, odd-parity p-wave pairing compatible with ferromagnetism is found in the vicinity of an orbital 
ordered phase. Furthermore, even-parity d-wave pairing with significant odd-frequency components is 
obtained. A possibility to observe such exotic superconductivity in manganites is also discussed briefly. 
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In the Bardeen-Cooper-Schrieffer theory for superconduc- 
tivity, 1 it was simply considered that magnetism suppresses 
superconductivity, since the singlet s-wave electron pair me- 
diated by phonon-induced attraction is easily destroyed by 
an applied magnetic field. However, since the pioneering dis- 
covery of superconductivity in Ce-based heavy-fermion ma- 
terial 2 and some uranium compounds, 3-7 it has been gradually 
recognized that anisotropic superconducting pair mediated by 
magnetic fluctuations generally appears in strongly correlated 
electron materials. In particular, due to successive discover- 
ies of superconductivity near an antiferromagnetic phase both 
in d- and /-electron systems, nowadays it is confirmed that 
"magnetism is a good friend to superconductivity". 8 

When we turn our attention to the relation between ferro- 
magnetism and superconductivity, it was discussed that criti- 
cal magnetic fluctuations can mediate triplet Cooper pair. 9 In 
fact, superconductivity has been observed in a ferromagnetic 
phase of uranium compounds such as UGe2, 10 URhGe, 11 
Ulr, 12 and UCoGe. 13 However, for /-electron systems, a mi- 
croscopic theory for superconductivity has not been satisfac- 
torily developed so far, mainly due to the difficulty in multi- 
orbital nature and strong spin-orbit coupling of / electrons. 

A way to overcome such a situation is to exploit a j-j 
coupling scheme. Along this research direction, the present 
author has developed microscopic /-electron theories on the 
basis of the j-j coupling scheme. 14-16 In the model, one /- 
electron state is characterized by an appropriate linear com- 
bination of the z component of total angular momentum j. 
Usually it is convenient to use the basis which diagonalizes 
the crystalline electric field (CEF) potential. In any case, we 
accommodate plural numbers of / electrons in such one-/- 
electron states due to the effect of Hund's rule interaction. 

For the case of cubic CEF potential, it is well known that 
the j=5/2 sextet is split into doublet and Tg quartet. Since 
the iy orbital has nodes along the cubic axes, it has strong lo- 
calized nature, while iy states have itinerant nature in com- 
parison with iy electrons. This orbital-dependent duality of 
/ electrons seems to be a key issue of rich phenomena in /- 
electron materials. Since electrons in localized iy and itiner- 
ant Tg orbitals are coupled with the Hund's rule interaction, 
we envisage a situation similar to double-exchange mangan- 
ites with mobile e g and localized t2 g electrons. 



In this Letter, a double-exchange scenario for the emer- 
gence of ferromagnetism in cubic uranium compounds is pro- 
posed on the basis of the orbital-dependent duality nature 
of 5/ electrons. We also propose some experiments to con- 
firm the double-exchange ferromagnetism in cubic uranium 
materials. In the ferromagnetic phase, we obtain the reduced 
Hamiltonian with active orbital degree of freedom. By analyz- 
ing the model within a random phase approximation (RPA), 
we find both odd-parity p-wave and even-parity d-wave pair- 
ing states in the vicinity of an orbital ordered state, suggest- 
ing orbital-related quantum critical phenomena. Finally, we 
briefly discuss a possibility of superconductivity in mangan- 
ites, which is well described by the double-exchange model. 

First we briefly explain the j-j coupling scheme. We in- 
clude the spin-orbit coupling so as to define the state labelled 
by the total angular momentum j, given by j=s+l, where 
s and £ are spin and angular momenta, respectively. For /- 
orbitals with £=3, we immediately obtain an octet with j=7/2 
and a sextet with j=5/2, which are well separated by the spin- 
orbit interaction. Since the octet level is higher than the sex- 
tet one, it is enough to consider j=5/2 sextet when local /- 
electron number is less than six. 

Next we define the one /-electron state in the cubic crystal 
structure. It is well known that under the cubic CEF potential, 
the sextet of j=5/2 is split into doublet and iy quartet. 
Note, however, that the ground state depends on the crystal 
structure. For instance, in the AuCu3-type cubic structure, the 
energy level for iy doublet is lower than that for Tg , while 
for CaF2-type cubic structure, Tg quartet becomes the ground 
state. In this paper, we assume the case with iy ground state. 

Since we consider the metallic uranium compounds, the 
valance of uranium ion takes the value between three and 
four, corresponding to the local /-electron number between 
three and two. When we accommodate two or three electrons 
in iy and Tg levels, we find two possibilities of low- and 
high- spin states, if we borrow the terminology of d-electron 
systems, depending on the balance between the Hund's rule 
interaction and the CEF splitting between and Tg levels. 

Readers may consider that the high- spin state is always sta- 
bilized in /-electron ions, but we should note that the effec- 
tive Hund's rule interaction J e ff in the j-j coupling scheme 
is reduced from the original Hund's rule coupling among /- 
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Fig. 1. Charge distribution of (a) T 8 a , (b) T 8 b , and (c) T 7 states. 



orbitals Jh as J e ^=Ju(gj — 1) 2 =Jh/49, 14 where gj is the 
Lande's g-factor and gj=6/7 for J=5/2. In fact, we have pro- 
posed the low-spin state for actinide ions to understand spin 
and orbital structure of AnTGas (An=U and Np; T=Ni, Pt, Fe 
and Co) 17 ' 18 and multipole order in NpC>2. 19,20 In this paper, 
on the other hand, we attempt to find new possibility of high- 
spin state concerning ferromagnetism and superconductivity. 

Now we discuss the /-electron kinetic term in a tight- 
binding approximation. When we evaluate /-electron hop- 
ping amplitude t* T , for nearest-neighbor hopping via the a 
bond between adjacent / orbitals, it is given by t% a =—V3t^ b 
=-V3t% a =3t&=3t/4, t y a =V3t y ab =Vst y ba =3t&=3t/4, and 
t* b =t, where indices a and b distinguishes two T 8 states (see 
Fig. 1) and t is given by t=3(f fa)/7 with the use of Slater- 
Koster integral (ffcr). 21,22 Note that orbital is localized, 
since the corresponding wavefunction has nodes along the 
axis directions, as shown in Fig. 1. On the other hand, Tg 
orbitals are itinerant and their hopping amplitudes are just the 
same as those of e g orbitals of 3d electrons, 15 ' 23 since is 
isomorphic to r3 x Tq, where Ts indicates E representation 
for the orbital part and Tq denotes the spin part. 

As mentioned above, we assume the high-spin state in this 
paper. Namely, the Hund's rule interaction works among 
and Tg orbitals. Note that Coulomb interaction in Tj states 
is larger than those for Tg ones in the order of J e ff. The 
difference of the magnitude of Coulomb interaction between 
itinerant and localized orbitals is not significant in compar- 
ison with d-electron systems, but in the combination with 
the orbital dependent duality nature, we arrive at the double- 
exchange model, which is used as a canonical model for man- 
ganites. 15,23 In this model, in order to gain the kinetic energy, 
the ferromagnetic phase appears, which is called the double- 
exchange ferromagnetism. This is established in the qualita- 
tive understanding of ferromagnetism in manganites. 

In order to confirm the emergence of double-exchange fer- 
romagnetism in uranium compounds, we propose a couple of 
experiments in analogy with manganites. One is the observa- 
tion of large negative magnetoresistance phenomenon. Since 
electrons can move smoothly in the ferromagnetic phase in 
comparison with the paramagnetic one, the resistivity is dras- 
tically decreased, when we apply a magnetic field on double- 
exchange materials. The large negative magnetoresistance in 
cubic uranium compounds may be an evidence for the double- 
exchange mechanism. For instance, in /3-US2, large magne- 
toresistance phenomenon has been observed, 24 although this 
material does not have cubic structure. 

Another is more direct evidence for the relation between 
the Curie temperature Tq and the kinetic energy. In mangan- 
ites with relatively wide bandwidth for conduction electrons, 
it has been observed that Tq is increased with the hole dop- 



ing. 23 Since the double-exchange ferromagnetism occurs so 
as to gain the kinetic energy, the ferromagnetic transition oc- 
curs more easily when electrons can move smoothly. Thus, we 
propose the appearance of the ferromagnetic metallic phase 
due to the application of hydrostatic pressure or the hole dop- 
ing on insulating and/or antiferromagnetic states of cubic ura- 
nium compounds. In the case of uranium compounds, hole 
doping can be done by thorium substitution. It is a drastic phe- 
nomenon that ferromagnetism appears due to thorium doping 
into antiferromagnetic uranium compounds. 

When Tg is lower than Tj with large Hund's rule cou- 
pling, the double-exchange ferromagnetism occurs for itin- 
erant Tj and localized Tg . Such a situation is realized in 
Nd-based filled skutterudite compounds, in which ferromag- 
netism is frequently observed. It may be interesting to seek for 
evidence of double-exchange ferromagnetism in such materi- 
als. Note also that the present mechanism cannot be directly 
applied to the tetragonal Uranium material, but it works even 
in the tetragonal system, when the CEF level splitting among 
Kramers doublets are less than J e ^ and the lower level has 
localized nature. 

Let us discuss the superconductivity in the ferromagnetic 
phase. For the purpose, we consider the spinless /-electron 
model with active orbital degree of freedom as 



H = 



i,a,T,T' 



(1) 



where fi T is the annihilation operator for an /-electron in the 
r-orbital of Tg at site i, m r =fl T fi r , and U is the inter-orbital 
Coulomb interaction. Throughout this paper, we set C/=4£, 
which is less than the bandwidth 6t. 

Note that if t" T ,=t5 TT r, H is equivalent to the well-known 
Hubbard model and we simply deduce that d-wave super- 
conductivity appears near the antiferro orbital-ordered phase. 
However, in actuality, electrons hop among different adjacent 
orbitals. The type of superconductivity in such a realistic mul- 
tiorbital system has been discussed actively, 25-27 and quite re- 
cently, it has attracted much attention due to the discovery of 
Fe-based superconductors. 28 

The non-interacting Green's function G is given by 



G-Hk) 



lUJ n 



M £-kaa 
~£kba 



l(x) n 



~£kab 

- H — Shbb 



(2) 



where we introduce the abbreviation fc=(fc, iuj n ), k is the mo- 
mentum, uo n =(2n + 1)ttT is the fermion Matsubara frequency 
with an integer n and a temperature T, Ehaa =3t(cosk x + 
cos k y )/2, Shbb =t(cos k x + cos k y + 4cos/c 2 )/2, e kab =e kba 
=— >/3t(cos k x — cos k y )/2, and a chemical potential fi con- 
trols the Tg electron number (n). Since Tj electron is as- 
sumed to be localized, the cases of (n)=l and 2 correspond to 
U 4+ and U 3+ ions, respectively. 

In order to discuss superconductivity, we solve the lin- 
earized gap equation for anomalous self-energy <j>, given by 

0t iT M = -TJ2 E K TlT2 ^(k,k , )^(k , ), (3) 

n' k' ,t[,T2 

where K(k, k') = V(k, k')G(k')G(-k') and V is given by 
V(k, k')=J+ Jx(k - k')[I - Jx(k - k')}- l J 



L-Lx{k + k')[I + L X {k + k')}- l L. 



(4) 
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Fig. 2. (Color online) Phase diagram for U =4t near the quantum critical 
region. Inset shows the whole phase diagram. 
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Fig. 3. (Color online) (a) Anomalous self-energies with odd-parity vs. n 
at fe=(7r/3, 7r/2, 37t/5). (b) Sign of 0i(fe,i7rT) with odd-parity on the 
Fermi surfaces. Red and blue colors denote plus and minus signs, respec- 
tively, (c) Sign of 4>3(k, iirT) with odd-parity on the Fermi surfaces. 



Here J a b,ab=JbaM =L aa,bb=Lbb,aa=U, I denotes unit ma- 
trix, and Xr 1 r 2 ,r3r 4 (<?) =~ Efc G T1T3 (k + q)G TAT2 (k). Here 
q=(q, v n ), q is the momentum, and v n =2nirT is the boson 
Matsubara frequency. In the calculation, we use a32x32x32 
lattice and 1024 Matsubara frequencies. 

In Fig. 2, we show the phase diagram in the (/i, T) plane. 
The boundary curve is determined from the divergence in 
the RPA susceptibility. The inset shows the whole phase di- 
agram: In the region I (0< u/t <0.85), the orbital ordered 
state appears. The ordering vector is Q=(tt, tt) at /i=0, 
but it is changed as (7r, 7r, S), where 6 is monotonically de- 
creased with the decrease of fi and it eventually becomes zero 
for ji/t >0.7. In the narrow region II (0.85< [ijt <1.05), 
we find Q=(5, 5, tt) with 5=11tt/16. In the region III (1.05< 
fi/t <1.76), Q=(S,S,S), where S=tt for T/t >0.1, while 
S < it fox T/t <0.1. 

On the analogy of anisotropic superconductivity near an 
antiferromagnetic critical point, we expect the appearance of 
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Fig. 4. (Color online) (a) Anomalous self-energies with even-parity vs. n 
at fe=(7r/3, 7r/2, 37t/5). (b) Sign of 4>i(k, iirT) with even-parity on the 
Fermi surfaces. Red and blue colors denote plus and minus signs, respec- 
tively, (c) Sign of 4>3(k, iirT) with even-parity on the Fermi surfaces. 



superconductivity when the orbital order is suppressed. In 
the present case, as shown in Fig. 2, there appears supercon- 
ducting pairing state due to orbital fluctuations with (5, 5, S) 
around at a quantum critical point fi/t ~ 1.76. Note that 
orbital is not the conserved quantity, since there exists non- 
zero hopping amplitude between different orbitals. Thus, it 
is meaningless to define orbital singlet and triplet by anal- 
ogy with spin singlet and triplet in the standard single-orbital 
Hubbard model. Here the superconducting pair is classified 
only by parity. In fact, we find that the superconducting state 
is labelled by even- and odd-parity, as shown in Fig. 2. 

We remark that even- and odd-frequency components 
are mixed in the present case. In order to understand 
this point, it is convenient to redefine the anomalous self- 
energy as 0l(fc)=0aa(*O, ^2(k)=(j) hh (k), <t>z(k) = [<j)ab{k) + 

</>ba(k)]/y/2, and 04(fe)=[0ab( fc ) - <j>ba{k)}/ 'V%. First we note 
that the relation (j)j{k)=—(j)j{—k) always holds for j=l~4, 
since it is due to the fermion property. The odd-parity 
solutions are characterized by <j)°(k,iuj n ) =—<j)°(—k,iuj n ) 
=(j)°(k,—iuj n ) for z=l~3 and (j)%(k,iuj n ) =—<p^(—k : iuj n ) 
=—(j)\{k, —iuo n ). Note that <j)\(k) has odd-frequency property. 
On the other hand, the even-parity solutions are character- 
ized by </>i(k, iuo n )= </>*(— k, iuo n )= -0f (ft, —iu n ) for z=l~3 
and 04 (ft, ioj n )= 04 (—ft, ioj n )= 0|(ft, —ioo n ). Note that 0f (&) 
with i=l~3 have odd-frequency properties. 

Let us first examine the odd-parity solution in the low- 
temperature region. In Fig. 3(a), we plot 0^'s vs. n of uo n . As 
mentioned above, 0i's for i=l~3 are even-frequency func- 
tions, while 04 is odd-frequency one. The absolute value of 
04 is relatively smaller than those of the even-frequency ones. 
Since 04 (A:) is exactly equal to the amplitude for the antisym- 
metric pair of electrons on different Fermi surfaces, the contri- 
bution of 04 is suppressed. In Figs. 3(b) and 3(c), we plot the 
signs of 0i (ft, inT) and 03 (ft, inT), respectively, on a couple 
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of Fermi surfaces in the first Brillouin zone. From the results, 
the pairing symmetry is found to be wave. We do not show 
the results for 02 and ^4, but 02 is similar to 0i and the mag- 
nitude of </>4 is small in comparison with other components. 
The node positions of 03 are different from those of 0i , but it 
is due to the difference in local symmetry of Tg a and Tg 6 . In 
fact, we find 03~0i (cos k x — cos k y ). 

In Fig. 4, we show the results for even-parity solution in 
the high-temperature region. The n dependence is depicted 
in Fig. 4(a). We find that 0^'s for i=l~3 are odd-frequency 
functions, while 04 is even-frequency one. Also in this case, 
the contribution of 04 is relatively small in comparison with 
other components. Thus, the even-parity solution is charac- 
terized by the odd-frequency components, leading to a way to 
observe peculiar odd-frequency pairing. 29 From Fig. 4(b), the 
gap function is found to be characterized by <i-wave. As ob- 
served in Fig. 4(c), 03 seems to be s-wave, due to the relation 
of 03 ~0i (cos k x — cos k y ) . 

Note that n is related to the valence of uranium ion. 
The critical point of corresponds to (n)^1.6, i.e., 

tj3.4+ jj ie w ^ t ] 1 f ^ superconducting region for the va- 
lence of uranium ion is the order of 0.01. Namely, the region 
is limited, but the value in the middle of U 3+ and U 4+ is re- 
alistic for actual uranium metallic compounds. Thus, we be- 
live that the superconductivity induced by orbital fluctuations 
could appear in ferromagnetic cubic uranium compounds. 

Here we mention a possibility to apply the theory to man- 
ganites, which are well described by the double-exchange 
model. The superconducting region corresponds to (n)«1.6, 
which denotes e g electron number for manganites. The sit- 
uation indicates 0.4 electrons per manganese ion form the 
particle-hole symmetry. Thus, the situation is close to the half- 
doped manganites with orbital ordering. In cubic manganites 
with relatively wide bandwidth, the metallic ferromagnetic 
phase is known to appear near the orbital ordering. The pat- 
tern of orbital ordering is different from the present one, but it 
is expected to observe superconductivity in manganites with 
high quality near half-doping. We note that the superconduct- 
ing transition temperature T c in Fig. 2 seems to be higher than 
that of the single-band Hubbard model within the same RPA. 
The stabilization of the even-parity solution due to signifi- 
cant odd-frequency components seems to be relevant to the 
increase T c . Since this point may open a new route to high-T c 
materials, further investigations will be required in future. 

Five comments are in order. (1) We have discussed orbital 
ordering and superconductivity in the ferromagnetic phase, 
but in order to confirm that the Curie temperature is higher 
than the orbital-ordering temperature and T c , it is necessary 
to estimate the magnitude of Coulomb interaction among / 
orbitals. This point is out of the scope of this paper, but it 
is one of future problems. (2) We have pointed out that iy 
becomes localized orbital when we take into account only a 
bond for / electron hopping. In general, hopping amplitudes 
through 7r and bonds appear and effective hoppings through 
ligand anions exist. Thus, Tf is not perfectly localized in ac- 
tual systems. However, we still believe that orbital dependent 
duality has an important starting point for the discussion on 
ferromagnetism and superconductivity. (3) We have ignored 
normal self-energy effects, but it is possible to include them, 
for instance, in the fluctuation-exchange (FLEX) approxima- 
tion. Without considering the vertex corrections, it overesti- 



mates the normal self-energy effect such as damping of quasi- 
particle, but in future, we can perform the FLEX calculation 
in the combination with dynamical mean-field approximation. 
(4) We have discussed superconductivity in the ferromagnetic 
phase from a microscopic viewpoint, but in actuality, it is nec- 
essary to consider how magnetic flux penetrates the system. If 
the magnetic flux forms some pattern such as the Abrikosov 
lattice, it indicates the ordering of localized iy electrons car- 
rying magnetic moments. This point may lead to an inter- 
esting possibility of the coupling between flux-lattice forma- 
tion and spin-orbital order. (5) We have proposed the spinless 
model, but from a realistic viwpoint, we should include both 
majority and minority spin bands. However, the minority spin 
band is virtually ignored, when minority spin density is so 
small that the intra-orbital Coulomb repulsion is effectively 
reduced in comparison with inter-orital Coulomb inetraction, 
indicating that orbital fluctuations dominate spin ones. 

In summary, we have proposed the double-exchange sce- 
nario for the emergence of ferromagnetism in cubic uranium 
compounds. We have found orbital-related quantum critical 
phenomena such as odd-parity p-wave and even-parity d- 
wave superconducting states in the vicinity of orbital-ordered 
phase. This orbital- fluctuation-induced superconductivity is 
expected to be found in ferromagnetic cubic uranium com- 
pounds and cubic perovskite manganites near the half-doping. 
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